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One reason for investigating Lightning Detection And Ranging (LDAR) is to validate data 
from the Optical Transient Detector (OTD). A Time-Of- Arrival (TOA) procedure may be used 
with radio wave portions of lightning signatures. An antenna network is in place at KSC [1 ]. 

Algorithms are available [2,3] which provide the advantages of considerably simplifying 
the numerical estimation of source location and a clear error analysis. These algorithms involve 
judicious differencing of measurement data which allows information retrieval by performing a 
linear inversion. Other algorithms [4,5] require numerically finding intersections of 
hyperboloids. Such algorithms are more involved computationally, and the error analysis is not 
straightforward. The algorithm in [3] provides for a clearly defined error analysis as given in [6], 

Efficiency in the retrieval of lightning location depends on the configuration of the 
antennas. Certain combinations of source locations and antenna configurations will result in the 
matrix of the linear system obtained using the algorithm in [3] having small eigenvalues. This 
will lead to a magnification of error. For example, a system of 4 antennas numbered from 1 to 4 
sequentially around the perimeter of a rectangle yields eigenvalues which are proportional to t n 
+ t M = t x - 1 2 + 1 3 - 1 4 , where tj is the time of arrival of a signal from the lightning event to the j th 
antenna. Letting Rj denote the distance from the source to the j th antenna and c be the speed of 
light, 

( 1 ) 


Curved transit paths due to refractive effects in the atmosphere are not considered. For brevity, 
cable time delays have been neglected. Source locations on either of the 2 planes which are 
perpendicular to the plane of the rectangle and bisect any line segment connecting 2 
consecutively numbered antennas yield t\i + ty 4 = 0. Such a configuration is “blind’’ to any 
sources on these planes. 

A better configuration places antennas at the centroid and vertices of an equilateral triangle. 
The antenna at the centroid will be designated by the index 1; others are numbered 2 through 4. 
This configuration has been referred to as “ideal” in [7], The eigenvalues for this configuration 
are proportional to t n + /j 3 + / 14 . This is equivalent to saying that the configuration is ineffective 
when the distance from the source to the antenna at the centroid is equal to the average of the 
distances from the source to the other 3 antennas. It turns out that this is never the case. The 
question of suitability of this particular arrangement of 4 coplanar antennas has been reduced to a 
problem in geometry. This report presents a solution. 

The desired result is obtained by seeking the extrema of the sum of distances from the 
vertices to a point on a hemisphere of radius R centered at the centroid of the triangle. Full 
generality is recovered by allowing R to be arbitrary. We set up a Cartesian coordinate system, 
with points specified by (x, y, z). Its origin is at the centroid. Vertices are located at 


D 2 =Dy, D 3 =-D 



( 2 ) 


the distance from the centroid to a vertex is D. For the present, assume R > D. It will be 
convenient to consider the distances from the vertices on an orthogonal coordinate system. 
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Vectors from antennas 1 through 4 to the source will be denoted by R, X, Y, Z respectively. Let 
0 denote the angle between D 2 and R. Employing the law of cosines, it may be shown that 


3 R 2 =X 2 + Y 2 + Z 2 - 3D 2 . 

(3) 

It will be necessary to extremize 


F(X,Y,Z) = X+Y + Z 

(4) 

subject to (3). 



Using Lagrange multipliers or any other suitable method, the only extremum in the first 

octant is X = Y = Z = R 2 + D 2 ; F = 3*Jr 2 + D 2 > 3 R. This turns out to be a local maximum. 
Minima must be found by looking along boundaries. These are not necessarily in any coordinate 
plane. All of X, Y, Z must be greater than or equal to R-D. Eliminating Z reduces the problem; 

f(X,Y) = X+Y+^3{R 2 +D 2 )-X 2 -y 2 (5) 

must be minimized requiring that none of X, Y, Z be negative. An additional constraint is 
obtained from (3), 

X 2 + Y 2 <l(R 2 +D 2 +RD). (6) 



This constraint does not sufficiently limit/; X and Y cannot both attain their minimum values of 
R-D at the same time. Boundaries are shown in Fig. 1. The smallest values of Y for a fixed X 
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occur when z = 0 (the lightning event is in the plane of the triangle). This may be demonstrated 
by writing 

Y 2 = R 2 + D 2 - 2R • D 3 . (7) 

Elevating the lightning event from the plane of the triangle while keeping X and R constant is 
equivalent to rotating X and R about D 2 . Here, D 3 has no z component and y is constant. The 
value of Y is minimized when x is as small as possible (at z = 0). 

The slope of a tangent to POQ is given by 

dY _ sin(2;r/3 - 6>) X 

dX~~ sin# Y ' 

It is never positive and goes from 0 at (V* 7 + D 2 +RD,R-Dj to — oo at 

- D, ^ R 2 + D 2 + RD). At the point O, X = Y, 6 = njl , and dYjdX = -1. The coordinates of 

this point are X - Y = yfR 2 + D 2 - ~RD. 

The gradient 

Vf(X, T) = (1 - XfZ)X + (1 - Y/Z) Y (9) 

points in the direction of greatest increase of f. The region in Fig. 1 can be thought of as a 
mountain with its summit at X = Y = 4~R~ 2 + D 2 (the point labeled S ). On the square region 

0 < X < V /? 2 + Id 2 , 0 < Y < yfR 2 + D 2 , we have Z > X and Z>Y. Except at S, the inequality is 
strict. The gradient points up and to the right. Below the line segment 
OS, rc/2) < 6 < 2zr/3, and X >Y. The opposite is true above OS. The mountain has a ridge along 
this line segment. Below and to the left of S, it is clear that as X and Y decrease, Z increases (the 
components of the gradient remain positive). Values of f beneath or to the left of POQ are 
smaller than values on this curve. 

It will not be necessary to investigate the values of/along PT or QT. These curves can be 
mapped onto POQ by a suitable relabeling of X, Y, Z. Symmetry may be further exploited to 
reduce the effort. An interchange of X and Y will map OP and OQ onto each other. 

It turns out to be difficult to parameterize Y in terms of X and obtain the minima along the 
bounding curve by differentiation. We will bound the values of f along OQ by those on a 
simpler polygonal boundary below OQ. 

We begin at Y= R- D and require 

f{X, R-D) = X + (R-D) + Z(X, R-D)>3R. (10) 

This leads to the inequalities 
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2X 2 - 2{2R + D)X + 2R 2 - D 2 + 2RD < 0, R - D<X<R + ^ + 1 D. (11) 

The largest value of X that needs to be considered is Jr 2 + D 2 + RD. Extrema along this line 
segment, denoted by QB (see Fig. 2), must be at the endpoints; 


/(Jr 2 +D 2 + /?£>, R- D^= R - D + 2 JY 2 + ~D 2 + RD > 3R, 
f(R-(yf3-\)D/2, R~d} = 3R. 


(12) 


The next portion of the bounding curve is drawn by leaving X = R-(j3 -\}D j2 and 

increasing Y to some point below OQ. The location of this point may be found by solving for 
cos# in X 2 = R 2 + D 1 - 2RDcos6. We find that 

Q V 3-1 Sd 

cos 6 = + , 
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(13) 


3 3/4 >/3 — 1 2 

< < — . 

2 3/2 4 3 

Using the law of cosines to write Y 2 in terms of R, D, and cos(2^/3 - 9), we see that Y decreases 
as cos(2;r/3 - #) increases. The smallest value of Y on POQ for X = R — {j3 — \^D j2 is 
bounded using (13); 

Y > Jit 2 +£> 2 - 47?D/3 > R- 2D/2 . (14) 

We can proceed up to point C whose coordinates are (r - (J3 -\^D j2, R -2D/Jj. 

As for the first part of the polygonal bounding curve, the extrema are at the endpoints. The 
minimum value off along BC is 3 R and occurs at B. The maximum is located at C and is given 
by 

2 R - [(V3 - l)/2 + 2/3]D + ^R 2 + (l4/9 + fityl) 2 + (V3 + 1/3 )RD > 3 R. (15) 

Proceeding as before, we fix Y=R- 2D/3 and see how far to the left we can go and still 
have / > 3 R. This time we obtain the inequalities 


V42/2-1 V42/2 + 1 ^ 

R -- — D<X<R + - — D. 


(16) 


3 3 

It will not be necessary to proceed from C all the way to the value indicated in (16). We can stop 
at a point E which has the same W-coordinate as O. From the estimate 
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f 


2 


(17) 


V42/2-1 

- 


<(R- 0.7 D) 2 < R 2 + D 2 -RD, 


we see that/X Y) > 3R at every point of this portion of the boundary. 

Finally, we proceed vertically from E to O. The derivative of / with respect to Y is positive 
everywhere along this line segment. The smallest value of f here is at E. 



Due to the relative positions of the actual and polygonal boundaries and the topography of 
the surface, the inequality is strict; X + Y + Z> 3R. When R > 0, the result for the case R < D 
may be inferred by interchanging R and D in the preceding calculations and showing X+ Y+ Z> 
3D > 3 R. The case R = 0 is immediate. 

The author would like to express his appreciation to W. Koshak, R. Blakeslee, H. Christian, 
and D. Phanord for their help and acknowledge the assistance of F. Six, G. Karr, L. Freeman, T. 
Shurtz, and S. Jefferson. 

REFERENCES 

1. Poehler, H. and C. Lennon, “Lightning detection and ranging system (LDAR), system 
description and performance objectives”, NASA Technical Memorandum 74105 (1979). 

2. Hager, W. and D. Wang, “An analysis of errors in the location, current, and velocity of 
lightning”, being reviewed for publication in J. Geophys. Res. (1995). 

3. Koshak, W., R. Solakiewicz, and H. Christian, “On the retrieval of lightning radio sources 
from time-of-arrival data”, in preparation for submission to J. Geophys. Res. (1995). 

4. Proctor, D., “A hyperbolic system for obtaining radio pictures of lightning”, J. Geophys Res. 
76(1971) 1478-1489. 

5. MacClement, W. and R. Murty, “VHF direction finder studies of lightning”, J. Appl. Meteor. 
17(1978)786-795. 

6. Twomey, S., Introduction to the Mathematics of Inversion in Remote Sensing and Indirect 
Measurements , Elsevier, New York, 1977. 

7. Poehler, H., “An accuracy analysis of the LDAR system”, NASA Contractor Report CR- 
154631 (1977). 


XLV-5 




